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ABSTRACT 

El 1 ipsometrlc measurements on plasma deposited "diamondl ike" amorphous 
carbon (a-C:H) films were taken in the visible, (E = 1.75 to 3.5 eV) The 
films were deposited on Si and their properties were varied using high tem- 
perature (up to 750 *C) anneals. The real (n) and imaginary (k> parts of 
the complex index of refraction N were obtained simultaneously. Follow- 
ing the theory of Forouhi and Bloomer (Phys. Rev. B34, 7018 (1986)), a 
least squares fit was used to find the dispersion relations n(E) and 
k(E) . Reasonably good fits were obtained, showing that the theory can be 
used'for a-C:H films. Morever, the value of the energy gap Eg obtained 
in this way was compared to the Eg. value using conventional Tauc plots 
and reasonably good agreement was obtained. 

INTRODUCTION 

The optical energy band gap of amorphous materials is usually found 
using a Tauc plot Cl], l.e., an extrapolation of (<xnE) 1/z versus E. Here 
a is the optical absorption coefficient, n is the refractive index and 
E Is the energy. In many cases, the refractive Index is almost constant 
or is unavailable, and (aE) 1/2 versus E is used (simplified Tauc). 

There are several drawbacks to the Tauc theory and procedure. First, the 
theory relates to absorption only, and cannot give the refractive index 
through the Kramers-Kroni g relation [2]. Second, there is ambiguity on. 
the energy range that the Tauc extrapolation is correct. Below a certain 
value of a, the absorption falls exponentially. This regime is called 
the Urbach edge [1]. The onset for this regime varies and values of order 
a = 10 3 cm- 1 [1] up to a = 10 4 cm- 1 [3] have been used. Empirically, many 
Tauc plots also start to deviate from a straight line at high a values 
[4,5] . 

Recently, an extension of the theory of optical absorption has been 
published for amorphous materials [2]. Later, the theory was extended to 
include more than one critical point, and was applied to crystalline semi- 
conductors [6]. In the Tauc plot derivation [1], the main assumptions are 
a constant matrix element and parabolic density of states for both the 
valence and the conduction bands. However, If the excited state has a 
finite lifetime x, the absorption probability has a damping factor [2] 
and Tauc plots are theoretically incorrect. The extinction coefficient k 
(k - <xch/2E where E is the energy) is calculated in [2,6] using the 
lifetime concept, obtaining. 


E - BE + C 

Here A is k(°°) and is proportional to M 2 /t where M is the position 
matrix element (M = < f j x 1 1 > ) ; B » 2 ( E c ^ cr 1 1 “ ^v.crit 1 ' where Eg.crit an<1 
E v crit are energies in the conduction and valence band respectively corre- 
sponding to a critical point, i.e., where k(E c>cr it - E v.crit> H-, 3 max1 " 
mum. C is related to the lifetime x through fi/x * (4C - BO 1 ' 3 and 
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Eg - E c bottom - E v,top- 1 - e - • the optical bandgap. Using the Kramers- 
Kronlg relation, the refractive index was obtained [2]: 


n(k) 


nO) 


8 o E 


E^ - BE + C 


( 2 ) 


Here B 0 and C 0 are related to A,B and C through simple algebraic 
formulas, and n(®) is a constant. 

In this paper we will examine if this new result applies to "diamond- 
like" carbon, also denoted a-C:H (amorphous hydrogenated carbon). If 
applicable, the result can give the n(E) function from a known absorption 
spectrum k(E) and a single refractive index measurement. Also, the func- 
tions n(E) and/or k(E) can be used when their analytical form is required, 
e.g., optimizing anti ref 1 ection filters. In addition, this paper will dis- 
cuss the meaning of the experimental constants, A,B, C, and Eg. 

Amorphous hydrogenated carbon material is made almost exclusively in 
form of thin films. The natural choice for an experimental technique is 
el 1 ipsometry . The multiple angle of incidence, multiple wavelengths (MAW) 
technique [7 to 9] was used. MAW gives n and k of the film simulta- 
neously at all wavelengths measured, without the use of either a known 
dispersion relation or application of the Kramers-Kronig analysis. This 
technique was commonly used to analyze various semiconductor multilayer 
structures [10,11] when n and k of the constituents were known. Here, 
the MAW technique is used in a spectroscopic way, to measure unknown n(E) 
and k(E). The smallest value of k that is accurately determined by 
ellipsometry is of order 0.005, corresponding to a = 1000 cm* 1 in the visi- 
ble. Thus, the Urbach edge regime is almost eliminated from this work. 
Reflection ellipsometry has the added benefit of measuring the thin film 
samples on any substrate, in contrast to transmission experiments. As the 
optical energy bandgap Eg depends on the substrate material [5], possibly 
due to differences in conductivity, this advantage of ellipsometry can be 
crucial for actual applications. 


EXPERIMENTAL 

Samples were prepared on 3 in. diameter Si substrates using a 30 kHz 
plasma deposition unit. The power P used was in the range 50 to 200 W. 
with a constant flow rate of 7xl0 -E m^/min. Details of the growth chamber 
are given in [12]. Several samples cut from the 150 W wafer were annealed 
in nitrogen for 10 sec using a rapid thermal anneal module [51. The 
rotating analyzer ellipsometer set-up [13] is essentially similar to the 
one described in [10]. At each angle of incidence, data was taken in the 
wavelength range 3500 to 7300 A, with 100 A intervals. Five angles of inci- 
dence were used, usually in the range 55“ to 75°. Each measurement con- 
sisted of 100 rev, with 72 points per rev, taken at a rate of -50 rev/sec. 
Background subtraction was done at each point. Calibration of the absolute 
value of the el 1 i psometri c parameters y and A was done using [14], The 
double grating monchromator was controlled by a IBM-AT personal computer 
that was also used for all data acquistion and analysis,, t and A were 
obtained by Fourier analysis. The MAW inversion process was done by 
minimizing experimental ^ and A (or tan y and cos A) versus their cal- 
culated values [91. The model used included only the substrate and a homo- 
geneous film. This is a reasonable approximation, as the interface of 
a-C:H on Si includes only the native oxide and -5 A SiC x [15], and our 
films were above 1000 A thick. 


RESULTS 

Seven samples were measured. Simplified and regular Tauc plots were 
done for all samples. A representative plot (50 W sample) is given in 
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Fig. 1, showing a slightly concave result versus the theory. As^the anneal- 
ing ' temperature and/or deposition power rises, the experimental Tauc func- 
tion first becomes a straight line and later turns convex. We kept the 
extrapolation range to the linear part of the experimental function. A sum- 
mary of the energy bandgaps obtained from regular and simplified Tauc^plots, 
together with sample thicknesses t, are given in Table I. Results of n(E) 
and k(E) for the 50 and 150 W, 600 °C samples are given in Figs 2 and 3 
respectively. Figure 2 shows a rising n versus E while Fig. 3 show a 
decreasing n(E). All other samples have an interim type of behavior. Grid 
least square fits to Eg. (1) were done to obtain A,B,C and Eg simulta- 
neously. With these 4 parameters fixed, the value of n(=») was set to get 
the best n<E) fits. The solid lines in Figs. 2 and 3 were calculated using 
these parameters. Results for all samples are given in Table II, while Eg 
is also displayed in Table I. 


TABLE I a-C : H OPTICAL ENERGY GAPS IN eV 


P, w 

isa 

1 

Sqr t ( aE ) 

Sqr t < d£n ) 

Grid 

50 


1490 

1 .90 

1.90 

1.64 

100 

■ I 

2305 

2.65 

2.65 

2.70 

150 

1 

2810 

2.36 

2.36 

2.23 

aoo 

■ 

2960 

2.02 

2.02 

1.65 

150 

■ 

2560 

2.05 

2.06 

1 .90 

150 

■ 

1320 

1 .08 

0.B9 

1 .03 

150 

■ 

950 

0.57 

0.39 

0.39 


TABLE II DISPERSION RELATION PARAMETERS FOR a-C:H FILMS 


P, w 

■ga 

A 

B, eV 

C , (eV)2 

Eg, eV 

n ( ® ) 

50 

KB 

.06 

6.39 

12.67 

1.64 

1 .63 

100 


.01 

6.46 

10.51 

2.70 

1 .70 

150 


.05 

5.90 

9.75 

2.23 

1 .67 

200 


.11 ' 

5.99 

11.73 

1 .65 

1 .69 

150 


.06 

5.10 

7.20 

1.90 

1.70 

150 

■ 

.40 

3.66 

4.22 

1.03 

1.58 

150 

m ^ 

.85 

1.21 

1.85 

0.39 

1 .71 
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FIGURE 3 * REFRACTIVE INDEX AND EXTINCTION COEFFICIENT 
OF THE 150 W, 600 °C SAMPLE. SOLID LINES ARE DUE TO 
EQUATIONS (1) AND (2). 


DISCUSSION 


Optical energy bandgaps Eg found using the two Tauc procedures are 
identical. However, Eg found using the dispersion relation fits are almost 
always smaller than the Tauc plots, although the differences are not very 
large. It Is impossible to rule which bandgap is the ''correct" one. The 
dispersion relation functions do take into account, at least qualitatively, 
the existence of absorption at lower energy than the regular Tauc regime, 
the Urbach tail. This can be seen by the very small slope of k(E) near 
Eg, as shown in Fig, 2. However, a(E) is expected to be purely exponential 
in the Urbach edge [1]. Thus, theoretically, the dispersion relations do 
not describe this regime. Table I also shows a decrease in Eg with 
increasing deposition power, (except the 50 W sample) and with increasing 
annealing temperature, in agreement with prior results [4,51. The thick- 
ness increases with power [16] and decreases with annealing temperature. 

The quality of the fits shown in Figs. 2 and 3 is quite good. The 
n(E) fits are the real test of the theory and the figures show a reasonable 
agreement between experiment and calculation. All other n(E) experimental 
curves show an equal or better fit to theory, as they include almost con- 
stant n(E) values. In addition, the k(E) fit is definitely better than 
the Tauc fits. 

The value of the A,8,C parameters, as shown in Table II, are essen- 
tially constant as function of deposition power, but they change signifi- 
cantly versus annealinq temperature. The value of the lifetime h/x, 
deduced from (4C - B 2 )' /2 , has quite an amount of scatter. However, all 
results are around h/x s 2 eV. This is a large value, but it is compara- 
ble with other amorphous materials [21. The values of A are in general 
lower than for other materials [21, denoting a smaller position matrix ele- 
ment M in a-C:H. There is an order of magnitude increase in M with 
heating to 750 g C. It would be interesting to correlate this increase with 
the changes in the composition and crystal 1 ivi ty of a-C:H. The value of 
B for the room temperature deposited samples, is equal to the lowest values 
obtained for other materials [2], and drops markedly with annealing tempera- 
ture. The value of B/2 is characteristic of the critical point bandgap 
[6], denoting a sharp decrease in this bandgap, in parallel with the sharp 
Eg drop. However, we did not see signs of a peak in k < E ) at E * B/2, as 
expected for a critical point. 
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CONCLUSION 


Dispersion relations suggested in [2] are obeyed, although the n(E> 
function does not have a perfect fit. Values of the parameters for the 
bandgaps 5 and Eg show reasonable agreement with Tauc plots and with 
prior results. The position matrix element increases with increasing 
annealing temperature. 
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